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ABSTRACT 



An analytical model of a jet injected normally from a flat 
plate into a uniform crossing flow was formulated to provide 
a simplified method of predicting the interference effects 
arising from the complex flow fields induced by ship bow 
thrusters. This model was an extension of previous work 
based upon a description of the jet as a series of distributed 
vortices. The analysis takes into account the position of 
the effective source of the jet and the blockage due to the 
presence of the jet in the crossflow. For representative 
jet-to-crossflow velocity ratios, the flow field and pressure 
distributions were calculated utilizing different combinations 
of effective source position and blockage. The accuracy of 
the model was evaluated by comparison with the available 
experimental data. Although good agreement was achieved for 
large portions of the interaction field, several regions were 
identified requiring further analytical description. 
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I. INTRODUCTION 



In recent years there has been an increase in the instal- 
lation and use of hull mounted, tunneled bow thrusters to 
improve ship maneuverability in restricted situations, such 
as dynamic position holding and confined waters docking. 

This increase has created renewed interest in the development 
of an ability to predict bow thruster performance. To date, 
there has been relatively little pre-installation design 
consideration regarding precise power requirements or ultimate 
operational performance. The limited testing that has been 
conducted /"~l-4 7* has resulted in a recurring major problem, 
namely, radical variations in effective turning moment with 
minor changes in ship speed. Chislett and Bjorheden / — 4_7 
noted that "an area of low pressure is created downstream 
of the discharging jet producing a resultant suction force 
that has a shifting center of action with changes in ship 
speed," thereby altering the effective moment produced by 
the thruster. The ability to predict the magnitude of this 
resultant suction force and its center of action would permit 
designers to accurately determine powering and control re- 
quirements for thrusters. 

This problem of radical variations in effective side 
forces and turning moments with only minor changes in ship 
speed may be better understood by considering the general 

* 

Numbers in brackets refer to references listed. 
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equations of motion of a ship moving in the horizontal 

plane of the sea. From /” 5_7 these equations are: 

m(u - rv) = X 

m(v - ru) = Y 

I r = N 
z 

where m = mass of the ship 

u = velocity in the x-direction 

v = velocity in the y-direction 

r = angular velocity 

I = mass moment of inertia 

X = excitation force (s) in x-direction 

Y = excitation force (s) in y-direction 

N = excitation torque (s) 

(dots indicate differentiation with respect to time) 

To minimize the complexity of the following brief analysis, 
all forces and torques will be referred to the shin's center 
of gravity. See Fig. 1. 

When these equations are transformed into their linearized 
version /” 5 T, they become 
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where the subscripts u f u, v, v, r, r f denote differentiation 
with respect to the variable subscripts, and u q is constant 
reference surge velocity. The left sides of these equations 
explicitly represent the velocity and acceleration dependent 
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reactions due to drag, added mass, added inertia and other 
coupled motions; while the right sides implicitly represent 
excitation force(s) and moment(s), which are in this analysis 
due to the thruster. Since this is a linearized formulation 
utilizing small angle perturbations, the excitation force in 
the x-direction is negligible compared to that in the 
y-direction : 

X << Y 

therefore X=0. The y-direction excitation force can be 
expressed as 

Y - Yj + Y i (1) 

where Yj = nuU^, the ideal jet thrust produced by the 
thruster, and Y^ is the induced force created by the inter- 
action of the thruster discharge with the flow past the ship. 
Utilizing experimental data from /” 4 _7 in the range, 

0.1 £ < 0.25, Y^ can be approximated by a linear function 

of ship speed 
Y . 

^ = 0.22 + k(§-) (2) 

j j 

where -3.0 < k < -1.5. Substituting a nominal value of 
k = -2.0 and Yj = m j U j into equation (2) gives 

Y i “ ^U^O.22 - 2.0 <§-)] (3) 

Substituting equation (3) into equation (1) results in 

Y = nij U j [ 1. 22 - 2.0(§-) ] 



11 



which explicitly couples the y-direction motion to the 
x-direction motion, normally uncoupled in the linear 
approximation. In addition to the coupling, in the relevant 
velocity ratio range, 0.1 < < 0.25, there is an appreciable 

reduction in the net thruster force due to motion in the 
x-direction. Furthermore, the moment produced by the thruster 
is a function of the net effective thrust, Y, and a character- 
istic moment arm, x fc : 

N t = x t Y 

As reported by Refs. /~"l-3 /, the degradation of net effective 
thrust is accompanied by a shift in the center of action of 
the induced force, Y^ , which results in a change in the char- 
acteristic moment arm. This shift in the thruster's net 
effective thrust center of action is a direct consequence of 
the coupling of x-direction motion with turning moment. 

This simplified analysis, which neglected all excitations 

other than the thruster, points out the complex, coupled mo- 

\ 

tions arising from the interaction of the thruster discharge 
and the flow past the ship. To fully understand this highly 
non-linear, three-dimensional flow situation extensive 
analytical and experimental work is required. This thesis 
will attempt to provide some understanding of the interaction 
of a jet discharging into a uniform cross-flow in order to 
develop some ability to predict the degradation of net effective 
thrust. 
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II. BACKGROUND 



A survey of the literature pertaining to the installation 
and operation of bow thrusters revealed a surprising lack of 
experimental data, considering the number of units currently 
installed. English /T,27 \ Stuntz and Taylor / — 3 and 
Chislett and Bjorheden /” 4 7 present some data from model tests 
conducted at various facilities, but in these works there is 
no formalized recording of full scale test or actual instal- 
lation test data. However, further investigation led to a 
related field: Vertical or Short Take-Off and Landing Aircraft 
(V/STOL) research. The concept of using air jets positioned 
in aircraft wings to provide vertical take-off and landing 
capability correlates extremely well with flow patterns pro- 
duced by a bow thruster. As reported by Bradbury and Wood 
/~ 6 _7 for incompressible flow, the jet (thruster discharge) 

path and induced flow are dependent mainly on the momentum 

2 2 

flux ratio, Pj U j /pU , and independent of Reynolds number. 
Margason /"~7_7 anci Gordier /~8 _7 also reported that the 
effective velocity ratio is the predominate characteristic 
in determining the path of the jet. Therefore, until more 
extensive ship or model data have been collected, the re- 
sults of V/STOL research are used as a basis for predictive 
theories. 

Although there has been a good deal of research conducted 
in the area of air jets in cross-flows, the majority of the 
work has been concerned with analytically or empirically 
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defining the jet trajectory and jet cross-section geometry. 
The interference effects of the interacting flows have re- 
ceived relatively little attention. A fluid jet injected 
into a crossing stream has previously been described as a 
turbulent, three-dimensional, highly non-linear flow problem 
which, even with considerable simplification, requires ex- 
tensive computer time to numerically solve the appropriate 
form of the Navier-Stokes equations /~ 9 _ 7 . Therefore, early 
investigators have sought simplified methods for determining 
the interaction effects, such as pressure and velocity dis- 
tributions, and have realized the necessity of knowing the 
geometry of the jet: trajectory and shape. 

Experimental data form the basis for most trajectory 
formulations. Either pressure or velocity measurements have 
been taken in the flow field encompassing the jet-crossflow 
interaction region and curves fit through points of maximum 
pressure or velocity. These curves were then compared with 
flow visualizations. In some instances, semi-empirical for- 
mulations were derived from the conservation laws combined 
with experimentally determined constants. Abramovich /”* 10 7 
qualitatively described the turbulent jet in a deflecting 
flow and presented some empirical methods for predicting 
the trajectory of the deflected jet. Jordinson /Tl 7 pre- 
sented trajectory data by recording and plotting contours 
of total pressure coefficients. Keffer and Baines /~12 _7 
presented experimental results from which it was determined 
that for various jet strengths, the jet trajectories could 
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be represented by a single function. Gordier /~ 3 _ 7 , the 
only reported source of water jet-water crossflow experiments, 
presented an empirical jet trajectory formulation based on 
a curve passed through experimental points of maximum total 
pressure. Sucec and Bowley /"~ 9_7 formulated an analytical 
expression for the jet trajectory utilizing previous experi- 
mental information and the assumption that the distributed 
pressure force and entrained momentum flux could be approxi- 
mated by an aerodynamic drag force. Margason / — 7_7 utilized 
flow visualization techniques and pressure measurements to 
formulate a trajectory equation. In addition, his paper also 
reviewed and compared the results of a number of other 
studies concerned with trajectory prediction. All the above 
predictive formulations produced trajectories within a range 
of uncertainty that can be attributed to expected experimental 
error, differences in test procedures and facilities, and in 
the case of analytical formulations, simplifying assumptions. 
In view of this, it was concluded that any one of them re- 
flected the current state of the art in trajectory prediction. 

The ability to predict the jet trajectory is of definite 
importance, but the primary objective of this jet-crossflow 
research is to provide insight into the changes, in the other- 
wise uniform flow, created by injecting a jet. Experiments 
conducted by Bradbury and Wood /~ 6 _ 7 , Vogler / — 13 _ 7 , Fearn 
and Weston / — 1 4 _ 7 , McMahon and Mosher / — 15 _ 7 » Kamotani and 
Greber / — 1 6_7 and others, provide measured values of the 
pressure distribution on the surface surrounding the jet 
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orifice which have been used to compare the accuracy of 
predictive models. During the past few years, a number of 
approximate methods to predict the pressure distribution has 
been developed: Wooler /T7 ,187 / Wu and Wright /19 7, Adler 

and Baron /” 20 7 and Schmitt / — 2^7* Again, because of the com- 
plex nature of the flow, all these models have utilized some 
degree of empiricism. Either the model is formulated using 
experimentally obtained trajectory expressions or the model 
is formulated using analytical expressions whose coefficients 
are selected to provide best fit with experimental data. 

The majority of the models are based on integral techniques 
with simplifying assumptions, such as (i) the representation 
of the entire flow as two-dimensional / 197 , (ii) the external 
flow is irrotaticna.l , uncompressible and steady-state /“ 20 7 or 
(iii) the flow is inviscid except that viscosity is the 
mechanism that leads to entrainment /T8 7. 

With the aid of these assumptions, the following approaches 
were made. Adler and Baron formulated their model by integrat- 
ing momentum equations to describe the jet mixing field, without 
using empirical trajectory data. However, numerous parameters 
were derived using empirical correlations from prior research. 
Their results provide satisfactory agreement with experiment, 
but from the quantity of empiricism involved, one might ex- 
pect such agreement. Wooler / — 1 7 7 followed by Wu and Wright 
/ 19_7 utilized blockage-sink representations where entrain- 
ment of crossflow fluid was handled analytically with suitably 
chosen coefficients. Again the results had satisfactory 
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agreement with experiment. However, an earlier effort by 
Wooler /~ 18 7 minimized empiricism to the extent that only 
an experimentally determined jet trajectory was necessary 
to complete the model. By thus restricting the number of 
empirical parameters, this model has considerable appeal 
from an engineering point of view. Since visualization of 
actual jet-crossflow situations has shown that the jet is 
deflected in the crossflow direction and forms two contra- 
rotating trailing vortices, this interference model was 
based on the representation of the jet by a distribution of 
vorticity. Incorporating arguments from the aerodynamic 
theory of lifting bodies this vorticity distribution was 
quantitatively described along the experimentally determined 
trajectory. By virtue of minimal empiricism and therefore 
greater appeal, this model was selected as the basis for 
continued investigation in developing a better predictive 
tool . 
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III. REVIEW OF WOOLER’S VORTEX MODEL 



Before any attempt to develop an improved model was 
made, a detailed re-derivation of the basis formulation was 
performed in order to: (1) gain a better understanding of 

the logic of formulation and (2) confirm the reported results, 
review of the distributed vortex method follows. 

One area of interest was the geometry of the problem. 

The empirical trajectory equation used was 

x/d = B ^cosh(z/Bd) - lj (4) 

where x is measured in the direction of the mainstream and 

z is measured in the direction of the exiting jet. The 

coefficient B was empirically determined from Jordinson /Tl7 

2 

and is equal to 0.19 (Uj/U) , where Uj is the jet velocity 
and U is the mainstream velocity. To further describe the 
problem) (see Fig. 2), a system of natural coordinates 
attached to the jet is adopted. Relative to the fixed 
coordinate system and an arbitrary point £x p ,y p j on t ^- e 
surface surrounding the jet orifice, the natural coordinates 
are given by 

K = - [zsing + (x-x p ) cos a ] 

n - Y p (5) 

C = zcosa - (x-x ) sina 

P 

where 5 is the coordinate tangent to the jet; C, the coordi- 
nate normal to the jet in the direction of the center of 



18 



curvature; n , the coordinate perpendicular to C and ? ; 
and a , the angle between the x and 5 directions. 

Another extremely important part of Wooler's formulation 
was the calculation of the distribution of vorticity within 
the jet. This distribution of vorticity was determined by 
adopting methods commonly found in the aerodynamic theory of 
lifting surfaces /” 227* According to this theory, regions of 
flow external to the jet are taken to be irrotational and the 
deflection of the jet is due to a purely inviscid mechanism. 
This inviscid mechanism is expressed as a balance between the 
forces due to pressure differences across an element of the 
jet and the centrifugal forces associated with the jet curva- 
ture. These pressure differences define corresponding velocity 
differences according to Bernoulli's Theorem. These resulting 
velocity differences, in turn, are related to circulation 
according to Kelvin's definition of circulation. Therefore, 
with the above assumptions, the circulation is related to 
jet curvature, and this relationship is 



where T is the circulation per unit length, ds, along the 
jet; R, the local radius of curvature of the jet; and d, the 
diameter of the jet orifice. Utilizing the empirical tra- 
jectory expression / — Eq. (4]_7 in the arc length derivative 
results in 




(6) 



ds _ _1 dz 

R Bd cosh (z/Bd) 



( 7 ) 
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Incorporating Eq. (7) into Eq. (6) and integrating over an 



element of jet extending from z^ to z^ results in the follow- 
ing expression for the circulation of a jet element 



K = r ud (^> 



-1 Z 2 -1 Z 1 
tan e /Bd - tan e /Bd 



J(8) 



In order to determine the flow field disturbance on the plate, 
the jet was divided into a number of vortex elements whose 
strengths were calculated using Eq. (8). These elements are 
the origins of characteristic horseshoe vortex systems where 
the cross member is bound in the jet and the trailing sides 
are, according to Wooler's assumption, tangent to the jet and 
separated by a distance, d, the diameter of the jet orifice 
(Fig. 3) . Errors associated with taking the trailing vortices 
tangent to the jet instead of containing them within the jet 
as in the actual flow situation were assumed by Wooler to have 
negligible effect on the velocity field along the plate. Each 
of these horseshoe vortex systems produces an interference 
velocity on the surface surrounding the jet orifice. To de- 
termine this interference velocity at any arbitrary point on 
the surface, the Law of Biot and Savart / — 2 3 7 




L 



is applied to each bound vortex and its associated trailing 

cl cl 

vortices. For the bound vortex, - /2 < L < /2, and for the 

trailing vortices, < L < « , were the intearation limits. 
These integrations produced closed-form solutions for the 
interference velocities due to individual jet elements. The 
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total interference velocity components due to the entire jet 
are determined as the sum of elemental contributions and 



are given by 




where N is the number of jet elements and and z^ + ^ are the 
endpoints of the general element. The parameters u, v and 
w are defined as geometric coefficients resulting from the 
integration along the jet. The subscripts 1, 2, 3 refer to 
the bound, left and right (as viewed from upstream) vortex 
contributions respectively. They are given by 



u, - 



1 C 2+ c 2 



a + h. 






( n + J s ) 2 -)/? 2 n 2 + (n- h ) 2 



v_ = 



2 *7/ TT 2 

C + (r)~ h) 



1 + 



<n- h ) 2 



v, - 



3 * ~7 



“ + 



x, + ( n + % ) ‘ -yj c 2+ c 2+ (n + 2 
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W, = 



l ■ 71 



v- + 



n + % 



- k 



? U 2 "^ 2 (n ♦ H ) 2 



? 2 + c 2 + ( n - *s> 2 



w = n - i 
*2 



C 2 + ( n- h ) 2 



1 + 



■yJT 2 + 5 2 + (n - %) 2 J 



n + % 

w 3 = 



C 2 + (n + h ) 2 



1+ v 



C 2 + c 2 + (n + *s ) 2 J 



and 5, C , and n are given by Eq. (5) for each point on the 
surface. After calculation of the interference velocity 
vector at a sufficient number of Doints on the surface, the 
pressure coefficinet at each point, C , is given by 



C P 



P - P 

0 

T 



i 2 , 2 

- l - q /u 



h pu 



where q 2 /U 2 = f P-t 2u 



(#)• 



The double u and v 



components are the result of the use of an image system to 
establish the surface as a solid boundary thereby negating 
the w component of velocity in this plane. Upon substitution 
and simplification, one has 



C 



P 




(9) 



After completion of the review and the incorporation 
of several necessary corrections to the published work /T8_7, 
the entire formulation was coded for computer solution. The 
resulting plots of pressure contours for representative 
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velocity ratios (Figs. 4,5,6) show fair agreement with the 
experiments of Bradbury and Wood / 6_7. However, these 
plots also point out serious deficiencies in the model. In 
the far field, in the arc 35° < 8 < 145° (6 being measured 
counter clockwise from the ray extending downstream from 
the jet origin) , the method represents the jet interference 
effects well, but in both the upstream and downstream 
portions of the surface, there is a definite lack of agree- 
ment. In the upstream area, the blockage effect due to the 
presence of the jet in the uniform flow is not predicted, 
while in the downstream regions, additional wake effects 
have not been taken into account. However, in view of the 
large area of good agreement with experiment, Wooler's 
distributed vorticity model was selected as the foundation 
for further investigation. It is apparent that this model 
gives a fair representation of the actual contra-rotating 
vortices associated with the jet-crossflow interaction. 
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IV. EXTENSION OF WOOLER MODEL 



Upon successful confirmation and correction of Wooler's 
distributed vortex model and as a step towards developing an 
improved model, it was felt necessary to verify the trailing 
vortex assumption, that these vortices are tangential to the 
jet rather than contained in the jet. Retaining the basic 
theoretical concepts for calculating the incremental circula- 
tion of a finite length of the jet, containment of the trailing 
vortices was approximated by successive conforming discretized 
vortex elements. These vortex elements were formed by restrict- 
ing the length of the trailing vortices to the linear distance 
between z^ and Z 2 / the arbitrarily-chosen end points of the 
finite length of jet, ds. In order to satisfy the Helmholtz 

vortex theorem, that a vortex cannot end in the fluid, the 

\ 

above horseshoe vortex system was closed by connecting the 
trailing arms with another bound, contra-rotating vortex. 

(See Fig. 7) By making this alteration, vortex loops of in- 
creasing incremental strength, as given by Eq. (8), can be 
made to conform to the curvature of the jet, as in the actual 
flow situation. This resulting vortex pattern is analogous 
to that created by a continuously accelerating wing /*"* 23 _7. 

As a wing accelerates, the strength of the produced vortices 
correspondingly increase, a situation very similar to the 
change in elemental vortex strength associated with the jet 
curvature . 
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Using this contained-vortex system, the interference 
velocity due to the jet at any arbitrary point on the plate 
was determined by applying the Biot-Savart Law (See 
Appendix A). In this application, it was necessary to 
approximate the trailing vortices as straight-line segments. 
However, as the sizes of the jet elements become small, the 
linear approximation can be made to conform to the curva- 
ture with an exactness that is only limited by numerical prac- 
ticalities. The total interference components due to the 
entire jet result from the summation of elemental contribu- 
tions and are given by 



u 
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1 

3 



•u. 



0 M 




i=l 
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-1 



k+1 + V 



z k+l z k 
e - e 



1 + e 



(z k+l + z k } ' 



where N is the number of jet elements and z^ and z^ + ^ are 
the endpoints of the general element. Again, the parameters 
u, v/ and w are defined as geometric coefficients resulting 
from applying the Biot-Savart Law around the element. The 
subscripts 1,2,3 and 4 refer to the forward and after bound 
and left and right trailing vortex contributions, respectively. 
These are given by 
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n - h 



l Vc 2 2 + Co 2 + (n + 55) 
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w. 
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? B2 


? B1 


c B 2 + (n - h) 2 


y«B2 + ? B 2 + <» - ’ s > 2 v 


' ? B1 + ? B 2 + (n - *=> 2 . 



w, 



n + h 


r 5 b 2 


? B1 


c B 2 + (n * is) 2 


y ? B2 + ? B 2 + <1 + ^l 2 -V 


'4l + ! B 2 - <" + l5)2 ] 



where the various £' s , c's and n ' s are given by Equations 
(Al, A2 , A3) of ADpendix A. With these interference velocity 
components at numerous points on the plate, the pressure 
coefficient, Eq. (9) , was determined for selected points 
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on the plate. The resulting contours of constant pressure 
coefficient are shown in Figs. 8, 9, and 10. 

These resulting plots show a somewhat different agree- 
ment with the experimental data /” 6 _7 than did the results 
of Wooler's model; however, the discrepancy can be directly 
attributed to the manner in which the trailing vortices are 
treated. In Wooler's model, the total interference velocity 
due to the trailing vortices results from the cumulative el- 
emental vortices acting at different distances relative to 
a point on the surface, while in the conforming vortex model 
the trailing vortex contribution results from the cumulative 
elemental vortex strength acting at a single radius, (See 
Fig. 11). Therefore, Wooler's assumption of treating the 

trailing vortices as being tangent to the trajectory intro- 
duced an error that under some conditions tended to improve 

agreement, where the conforming vortex model removed this 
error at the expense of some loss in agreement, at least at 
the larger values of Uj/U. However, there have been 
studies /"” 12,217 that present the concept of an effective 
origin of the jet vortex system - a region where the turbu- 
lent mixing processes extend across the entire jet causing 
definite deflection of the jet and establishment of the contra- 
rotating vortices. Schmidt / — 21 7 formulated an empirical 
expression for the position of the effective source as a 
function of the jet-to-mainstream velocity ratio given by 



where a = U/U^ 



Z = (15 o/a) ~ 1 
e 

Applying this effective source expression 
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to the conforming vortex model produces better agreement 
with experiment (Fig. 12) . 

Comparison of the conforming vortex model, including an 
effective source, with experimental data shows the same, 
if not better, far-field agreement with an improvement in 
the predicted pressure distribution in upstream and downstream 
regions. Therefore, this model which more closely represents 
the actual flow situation, replaced Wooler's model as the 
basis for further investigation and improvement. 
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V. IMPROVEMENT OF CONFORMING VORTEX MODEL 



Evaluation of the conforming vortex model results in- 
dicated the continued lack of agreement with experiment in 
two distinct areas: upstream and wake. In order to improve 

the agreement in these areas, it was felt that additional 
corrections could be formulated from a qualitative analysis 
of the flow phenomenon. The discrepancies in the upstream 
area can be, at least partially, attributed to the blockage 
of the mainstream flow by the presence of the jet, while the 
discrepancies in the wake area are caused by turbulent flow 
separation, vortex shedding and turbulent entrainment. In 
view of the present uncertainty concerning the exact character 
of the wake area, it was decided that the development of 
corrections for the wake must await future investigation and 
experimentation. Therefore, the upstream blockage correction 
was undertaken. 

Near the point of injection, the physical appearance of 
a circular jet injected into a crossing flow is that of a 
cylinder. Interpreting the jet as a cylinder in an irrotational 
flow, potential flow theory predicts perturbations to the 
uniform flow that are given by 



u 

U 




cos 2 0 



v 

U 




sin 20 



( 10 ) 
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where r and 9 are cylindrical coordinates of points in the 
flow field and a is the radius of the cylinder. These 
velocity perturbations caused by blockage can be combined 
with the interference velocity components derived from the 
vorticity formulation to adjust the coefficient of pressure 
(Eq. 9) at arbitrary points on the surface. The initial 
(and most simple) assumption was to make the diameter of the 
blockage cylinder the same as the diameter of the jet orifice. 
The resulting pressure distribution indicated an excessive 
blockage that is not present in the real flow. The actual 
flow does not produce a solid interface, as would a cylinder, 
but instead an entraining interface surrounding a core of 
uniform velocity approximately equal to U ^ , the jet exit 
velocity. To better describe this, a potential cylinder of 
diameter less than the diameter of the jet was introduced. 
Using many combinations of cylinder diameters with different 
jet-to-mainstream velocity ratios, an expression for a in 
Eq. (10) was formulated in terms of the velocity ratio. 

a = 0.96 So 

Applying this formulation to the blockage-cylinder perturba- 
tion velocity calculations results in pressure contour plots 
(Fig. 13,14,15) that show improved agreement with experiment 
in the far field and upstream areas for jet-to-mainstream 
velocity ratios, 4 < -i < 8. However, as the velocity ratio 
increases, the model results and experiment diverge. The 
apparent cause of this disagreement can be resolved by 
qualitatively examining the actual flow. With increasing 
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velocity ratios, the entrainment of crossflow fluid by the 
jet becomes the more dominant interference factor, while the 
blockage factor is relatively less important. Wu and Wright 
/ 19_7 have concluded that "the amount of crosswind fluid en- 
trained increases with increasing jet speed" when the cross- 
wind speed is held constant. In addition, this increased 
entrainment appears to be especially important in the down- 
steam wake region and causes a further reduction in the 
blockage effect as compared to that due to a solid cylinder 
in a crossflow. The resulting effect of the increased rate 
of entrainment at higher jet speeds is larger interference 
velocity perturbations in the flow field surrounding the jet 
orifice near the plate. In view of this, the lack of agree- 
ment at higher velocity ratios can be explained while also 
pointing out the need for additional correction (s ) to account 
for the change in the relative importance of blockage and 
entrainment with changes in velocity ratio. 
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VI. CURRENT MODEL CAPABILITY 



Realizing the limitations of the current analytical 
model, it is nevertheless appropriate to investigate its 
predictive capability. Figure 16 shows the results of the 
analysis expressed as the force induced on a flat plate due 
to a jet injected at right angles to a crossflow. The force 
is expressed as a fraction of the ideal jet thrust and the 
negative values indicate suction forces (in a direction op- 
posite to the ideal jet thrust) . Correlation with experiment 
/“4_7 shows a difference in the value of k in Eq. (2). The 
value of k - -4.4 resulting from this analysis, although not 
in the predicted range, -3.0 < k < -1.5, still displays the 
correct trend and provides encouragement as to the potential 
success of a model based upon the present methods. The main 
difference in the approximation constant results from the 
lack of proper accounting for the entire entrainment effect, 
a factor which progressively degrades the model at higher 
velocity ratios. However, as an initial design approximation 
for low velocity bow thrusters active at relatively high 
ship speeds, the conforming vortex model presented here pro- 
vides a satisfactory prediction of the interference effects 
— a step toward improving the accuracy of oowering and con- 
trol determinations. 
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VII. CONCLUSIONS 



The conforming vortex model, modified with corrections 
based on experimental observations, is presented as an 
initial formulation of a simple yet physically consistent 
representation of jet-crossflow interference effects. Pres- 
sure distributions in the far-field are adequately predicted 
and the calculated results obtained with the model provide 
reasonably good agreement with experiment for the lower 
jet-to-crossflow velocity ratios. Agreement is less than 
satisfactory for the higher velocity ratios. However, the 
model does indicate the relative importance of different 
flow phenomena with changing velocity ratios. The most evi- 
dent weakness of the model is the lack of an adequate method 
for including entrainment and blockage effects near to the 
origin of the jet at higher velocity ratios. By restricting 
the range of application, the conforming vortex model demon- 
strates sufficient agreement to be utilized as a first 
approximation technique for designing bow thrusters. 
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VIII. RECOMMENDATIONS 



Further experiments should be conducted to provide more 
correlation data, especially in the case of water jet into 
water crossflow. Empirical constants utilized in the current 
model and those that will become necessary when trying to 
model increased entrainment can be refined when a larger data 
base is available. The encouraging results of this method 
of analyzing the flow field should not be overlooked, and 
additional corrections, as indicated in this analysis, should 
be made to improve its accuracy. 
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Figure 1. Sketch of bow thruster-crossflow interaction. 
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Figure 2. Description of coordinate systems 




Figure 3. 



Wooler Horseshoe Vortex System 
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Figure 4. Comparison Of Pressure Coefficient Contours 

Experiment / 6 7 , Wooler 

Velocity Ratio = 4 
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Figure 5. Comparison Of Pressure Coefficient Contours 

Experiment / — 6_7* Wooler — 

Uj/U = 8 
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Figure 6. 



Comparison 
Experiment 
Uj/U = 11.3 



Of Pressure Coefficient Contours 
/ 6 7 , Wooler 
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Figure 7. Conforming Vortex Loop System 
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Figure 8. Comparison Of Pressure Coefficient Contours 
Experiment ---/ 6 ~ 7 , Conforming Vortex Model 
Uj/U = 4.0. “ ~ 
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Figure 9. Comparison Of Pressure Coefficient Contours 
Experiment / 6 _ 7 , Conforming Vortex Model 

Uj/U = 8.0 
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Figure 10. 



Comparison Of Pressure Coefficient Contours 

Experiment / 6_7, Conforming Vortex Model 

Uj/U = 11.3 
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Figure 11. Comparison Of Wooler (a) And Conforming Vortex 

Model (b) Effective Radii For Induced Velocities 
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Figure 12. Comparison Of Pressure Coefficient Contours 

Experiment / 6_7, Conforming Vortex Model 

With Effective Source Correction .Uj/U =8.0 
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Figure 13. Comparison Of Pressure Coefficient Contours 

Experiment / 6 , Fully Corrected Conforming 

Vortex Model — .U./u = 4.0 
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Figure 14. Comparison Of Pressure Coefficient Contours 

Experiment — / 6 7» Fully Corrected Conforming 
Vortex Model — . U.7u = 8.0 
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Figure 15. Comparison Of Pressure Coefficient Contours 

Experiment / 6 7 , Fully Corrected Conforming 

Vortex Model .0 .7u = 11.3 
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Figure 16. Variation In The Induced Force Due To Jet-Crossflow 
Interaction 
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APPENDIX A 



FORMULATION OF INTERFERENCE VELOCITY COMPONENTS 



The formulation of the interference velocity due to the 
vortex loop system of the conforming vortex model is given 
below. The formulation requires the use of three separate 
coordinate systems that are functions of the endpoints of a 
jet element and the position of an arbitrary point on the 
plate. From Fig. (Al) , these are given by: 



= - ^z 1 sina 1 + (x 1 - x p ) cosa^j 



n = y r 



c = z.cosa, - (x. - x ) sina. 

]_ 1 1 1 p 1 



for the upstream bound vortex 



^2 = " [ Z 2 sina 2 + (x 2 " X p } COSa 2] 



n = y r 



= Z2Cosa2 - (X2 - x^) sini^ 



for the downstream bound vortex, and 



hi = ■ [ Z i sina B + (x i ’ V COSa B] 



n = y 

C B = z iC osa B - (x ± - x n ) sina B 



(Al) 



(A2 ) 



(A3) 
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for the trailing vortices, where x and y are the 

P P 

coordinates of an arbitrary point on the plate and x^, z^ 
X2»Z2 are the coordinates of the endpoints of an element 
of jet. The angles a^, and ct B are the angles between the 
x-direction and respective £'s and are given by 

= tan ^ ^csch (z^/B)j 
= tan ^ j^csch (Z2/B)j 
a B = tan -1 ^(z 2 - z 1 )/(x 2 - x 1 )j 



To determine the induced velocity at an arbitrary point, 
the law of Biot and Savart 



T f sing 

4 W7" 



L 



as 



(A4) 



must be individually applied to both the bound and trailing 
vorticies. Due to similarities in the integrations, only 
one bound and one trailing vortex formulation will be pre- 
sented and necessary changes applied for the other vortices. 
The circulation of a finite amount of jet, given by /~~ Eq. 



( 8)7 



K = 



7T 

2 



Zi 

U 



,v 


-1 Z 2 /Bd 


, z./Bd 


j 


tan e 


- tan e 



is the incremental contribution which is added to each 
successive element, such that the total circulation for 
a given element is given by 
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n 




■ -i z k+i /Bd „ -i z k /Bd 
tan e - tan e 



(A5) 



Therefore, the induced velocity due to any of the vortex 
filaments associated with the n-th element, 1 £ n N, 
where N is the number of jet elements, can be determined 
from 



I ^ f 

! n' ~ 17 J 



sinB 



ds 



(A6) 



where L is the length of the filament, B is the angle 
between the vortex filament and the radius, r, to an 
arbitrary point on the plate. 

For the forward bound vortex (See Fig. (A2)) 



sinB = 



V^ 2 + 



V?1 2 + C x 2 + (n + 



(n + s) 



and 



2 2 2 
= + (rj + s) 



Substituting these equations into Eq. (A6) gives 

r 



l nl' 



n 

17 



2 2 
+ ? i 





[c x 2 + ? x 2 + (n + S) 2 ] 



T\ 3/2 



ds 



(A7) 



the magnitude of the induced velocity at an arbitrary point 

on the plate due to the bound vortex of the n-th element. 
The x-direction component of this velocity is given by 
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q xl = I q nl I [cosy 1 cosa 1 - siny^sinaj 



(A8) 



where cosy. = 



Vq 2 + C-l 2 



and siny, = 






Inserting Eq. (A7) into Eq. (A8) and applying the integration 
limits gives 

r Ji 

r 

rr n 

q xl 4 it 

1^5 



2 2 2 
+ (n + s) z 



cos ct^ - 



h 



(A9) 



2 2 2 

5i + Ci + (n + s) z 



37 T“ ds sinct i 



Integration of Eq. (A9) results in 



‘xl 



n 

4 TT 



n + h 



~~2 I 

h + 'i 



- h 






2 2 
+ + 



(n + %) 



V? 



+ + (n-*s) 






qosa. 




n + h 



T 

’1 



+ 57 + (n + %) 






+ % c 

+ C 2 + (n-%) 2 



:ina. 



(A10) 



the x-direction component of induced velocity due to the 
forward bound vortex. Similarly, the x-direction component 
of induced velocity due to the downstream bound vortex is 
given by 
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n 



*x2 TTt 



n + H 



n - h \cosg. 



- ? 2 + c 2 V?2 2 + c 2 2 + (n + * 5 ) 2 V? 2 2 +C 2 2 +(n-! 5 ) 2 y 



n + '*5 



n - *s 



«2 + «2V 2 S * c 2 2+ <" + « 2 ' V« 2 2 *c 2 2 +<n-«?' 



sina. 



(All) 



The approximation of the bound vortices as being linear 
eliminates the y-direction component of induced velocity 
and since an image system is utilized to formulate the flat 
plate, it is unnecessary to calculate any z-direction 
component of induced velocity. » 

For the left trailing vortex filament, as viewed from 
upstream (See Fig. (A3) ) . 

Vc B 2 ♦ (n - *) 2 ' 

sinB = ======_ 

Vc B 2 + ( n- h) 2 + s 2 

r 2 = Cg 2 + (n - h) 2 + s 2 



Inserting these expressions into Eq. (A6) gives 

r 



*n3 



n 
4 IT 



C B 2 + (n - h ) 2 



k 2 



+ ( n - h)* + S 



1 



371 



ds 



(A12 ) 



the magnitude of the induced velocity at an arbitrary 
point on the plate due to the left trailing vortex of the 
n-th element. Resolving this velocity into x and y com- 
ponents gives 
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l x3 



and 



y 3 



J n3 



cosy D sina_ 

O O 



(A13) 



*n3 



siny 



B 



where cosy B = 



n - h 



C B 2 + (n - h) 2 



and siny B = 



’B 



C B 2 + (n - h) 2 



Inserting Eq. (A12) into Eq. (A13) and applying integration 
limits gives 

r y O 



x3 



F n f B2 

4? / r 2 . H, . 2 ~ 213/2 ds Slna B 

^ B 1 B J 



y3 



F n f B2 

= ^ J, [PP 



’B 






•Bl l B 

Integration of Eq. (A14) yields 



ds 



(A14) 



_ _ n 
*x3 4 tt 



n - h 



’32 



1 7 , , 

l ? b + (n - h) VL 2 



’Bl 



sina 



, ? B ‘+(n-> i ) 2 +C^ 2 V; B 2 +(n->» 2 t£ 2 i 



3 



n 



^y3 ^ 4 tt 



’B 



’B2 



’Bl 



L«B + ^“» VVc B 2 + (n-« 2 + C 2 2 Vc b 2+ (n— % ) 2 +£ B1 /J 



(A15) 

the x and y corroonents of induced velocity due to the left 
trailing vortex. Similarly the x and y components of the 
induced velocity due to the right trailing vortex are given 
by 
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n 



q x4 4 tt 



n + h 



’B2 



’B1 
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l ? b + (n + *s) \Vc B 2 +(n+^) 2 +5 B2 Vc b 2 + (n+%) 2 +? B1 / 



sina T 



n 



*y 4 4 77 



>B2 



*Bl 



c b + (n + %) Vc b 2 +( n+%) 2 +C B2 Vc b 2 + (n+5s) 2 +? 2 i/-J 



(A16 ) 

Therefore, the total x and y components of induced velocity 
due to an entire element of jet are 



0*T ’ Oxl + q x2 + 0 x3 + 1 xi 
*3yT “ q y3 + q y4 



(A17) 



Inserting Eqs. (A10) , (All) , (A15) , and (A16) into Eq. (A14) 

yields 

r n r 1 

q xT = 4 ^jf ! cosci^ - u 2 cosa 2 -w 1 sina 1 +w 2 sina 2 + («w 3 +w 4 ) sinctg | 



q yT = Tff [ V 3 " V 4] 



(A18) 



where 
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n + h 



n - h 
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w, = 



1 - 2 . : : 

*1 + c i 



n + h 



n - h 



w 2 T~1 



c 2 + ? 2 - 



-V?i 2 +Ci 2 + ( n+>i ) 2 Ve 1 2 +C 1 2 + <n -%) 2 -* 

n + % n - % 
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S B + (n-*s) 



’B2 
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•Vd 0 +c 0 2 +(r l-is ) 2 Vc Bl +s B 2 +( n-Js) 



w 4 ~ 2 



n + *5 



c B + (n+h) 
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Bl S B 



Substituting the expression for T n from Eq. (A5) into Eqs. 
(A18) and non-dimensionalizing yields 
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U 



l i 

8 V U 




-w^sina^+w 2 sina_^ + ^+ (-w^+v^) sin ctg+u^cos 
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k=l 
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,1 + e 
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Figure Al. Coordinate Systems For Conforming Vortex Model 
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Figure A2. Geometry Of Bound Vortex Integration 
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Figure A3. Geometry Of Trailing Vortex Integration 
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E = C .1<3*VR**2 

CALCULATION OF INTERFERENCE VELOCITY 
L = C 

CALCON = DFLOAT(NANG) - 1.0 

C INCREMENT RADIUS (INFLT VALUES) 

CC 1 AO I = 1 , N RAD 
RCALC = R (I ) 

C INCREMENT THETA (5 DEGREE INCREMENTS) 

CO 120 J = 1 i NANG 
K = J-l 

THETA ( J ) = DFlQAT(K)* (PI/CALCON) 

C SHIFT TO CARTESIAN CCCPC1NATES 
L = L +1 

COST = DCOS ( T HETA ( J ) ) 

SINT = DS IN (THETA( J ) ) 

XP ( L ) = RCA LC *COST 
YP ( L ) » R CALC ’♦SI NT 
C 

21 = 0.0 
22 = C.00CC1 

HKSAV = 0.0 

ICC 21 = 22 

22 =» Z1 + ZINC 
211 = Z1 ♦ CELTAZ 
111 ~ 12 * DELTAZ 

>1 = 8*(DC0SH(Z1/B)-1.0 ) 

>2 * e*(DCOSH(Z2/B)-1.0) 

TANAe = ( 12-1 1 ) /( X2-X1 ) 

TAN A 1 = 1 .0/C S INH ( Z 1/8 ) 

TANA2 = 1 .0/0 S INH ( 22/ E ) 

ALFCNE = CATAN(TANA1 ) 

ALFT V*0 = CATAMTANA2) 

ALFEAR = DATAN (TANAE ) 

SINAI = DSI N ( ALFONE ) 

SINA2 = DSIN(ALFTWC) 

SINAE = CSIN( ALFBAR ) 

C0SA1 = CCCS( ALFONE ) 

CCSA2 = OCOS ( ALFTWO ) 

COSAe = DCGS ( ALFBAR ) 

C NATURAL COORDINATES IN TERMS OF CARTESIAN COORDINATES 

XICNE = -(Z1HSINA1 + ( Xl-XP(L) )*CCSA1) 

XITVC = -(Z22+SINA2 + ( X2-XP( L ) ) *CCS A2 ) 

XIBAR1 = -( 211*SINAE + (Xl-XP (L ) ) *CCSAB ) 

XIEAR2 = -<222*SINAE + (X2-XP (L ) ) *CCSAB) 

ETA = YP ( L ) 

2ETA1 = <Z11*CCSA1 - ( X 1-XP ( L ) ) *S IN A 1 ) 

ZETA2 = ( Z22* COS A2 - ( X 2-XP ( L ) ) *S IN A 2 ) 

2ETAe = (211 + CCSAB - (Xl-XP (L ) ) *S INAE ) 
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C GEOMETRIC CONSTA MS 

ETPl = ETA * C.5 
ETMI = ETA - C.5 
SA = XICNE**2 
S £ = XITWG**2 
SC = X I B A R 1 ** 2 
SC = XIBAR2**2 
SE = ZETA 1**2 
SF = ZETA 2**2 
SG = ZET A E**2 
SF = ETPL**2 
SI = ETMI **2 

CA = DS CRT ( SA + SE + SH) 

CB = DSQR T < SA + SE -» SI) 

CC = OSCRT(Se + SF 4 SH) 

CE = DSQR T( SB 4 SF 4 SI ) 

CF = DSQR T ( SD 4 SG 4 SI ) 

£G = CS CRT (SC + SG 4 SI ) 

CH = DSQRT (SC + SG 4 SH ) 

Cl = CS CRT ( SC 4 SG 4 SH) 

CCM = (ETPL/CA - ETHI/CB)/(SA 4 SE) 

CCN2 = (ETPL/CC - ETNI/CE)/(SE 4 <F) 

CGN3 = (XI8AR2/0F - X I B AR 1/CG ) / ( S G 4 SI) 
CCN4 = ( X I8AR2/0H - XIBAR1/DI )/(SG 4 SH) 

LI = ZETA 1*CC M 

L2 = ZETA2*C0N2 

V 2 = ZETA6*CCN3 
V4 = ZETA B*CO N4 
VI = XICNE*COM 

V 2 = XITWC*C0N2 
V3 = ETMI *C CN !'• 

V4 = ET PL*CCN 4 



C INCREMENTAL VORTICITY TERM 

FF = <(DEXP(Z2/B)-DEXP(Z1/B))/(1.C4CEXP((Z24Z1)/E)) ) 

c 

FK = <1. 0/8.0)* VR**2 * OATAN(FF) 

FK = HK 4 HKSAV 
FKSAV = HK 

C INCREMENTAL INTERFERENCE VELOCITY COMPONENTS 
LP * C.O 

LP = FK*( ( — Vtl4$INAl4V(2*SINA2+(“Vl3+ViA )*SINA8 + L1*CCSA1 ) 
X -(L2*COSA2)J 
C 

VP = 0.0 

VP = HK* { V3 - V4) 

C 

C 

U ItJ) = U(I* J) + UP 
V(I ,J) = v< ItJ) + VP 

C TEST FCR CONVERGENCE 

TEST = CABS(UP) 4 DABS(VP) 

IF (TEST. LE.C.C000C1 ) GC TO 110 
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c 

GC 1C 100 

c 

110 CONTINUE 

C CALCULATION OF BLOCKAGE VELOCITY TERN 

IF(NTYPE.EQ.C) GO TC 120 

IF < NTYPE. EC .2 ) GO TC 120 
C 

R 6 LC C K = C.96/CSQPT (VR ) 

CB LOCK = 2.0*PBL0CK*C 

PK = <<RBLCCK/R< I) )**2)/2.0 
LC = RK*(SINTx*2 - CCST**2) 

VC = -2 • 0*RK* S INT*CCST 
UIiJ) = U( I, J) + UC 
V ( I ,J) = V( I, Jl + VC 

C CALCULATION CF PRESSURE COEFFICIENT 

120 CP<I,J) = -4 .C*( U ( I * J ) + U(I»J)**2 ♦ V < I » J ) * *2 ) 

C CHANGE CES IGNAT I CN OF VARIABLE 

C 

Z P < L ) = CF( I« J) 

C 

120 CONTINUE 
140 CONTINUE 

C CALCULATION OF INCREMENTAL FORCE 

RO = 0.0 
P I = 0.0 
AREA * 0.0 
PAREA = 0.0 
FCPCE = 0.0 
YMON = 0.0 
TCTFCR = 0.0 
TYMCM = 0.0 

NAPEA = NFAC - 1 

CC 1 7 C I = lt NA FEA 
K = I + 1 

pi = R<n 
RG = R ( K I 

C CALCULATION CF INCREMENTAL AREA 



AREA = PI *( R0**2 - P 1**2 ) 

FAR EA = AREA/ <2.0*CALC0N) 

CO 160 J=1,NANG 

C CALCULATION CF AVERAGE PRESSURE COEFFICIENT 



CFAVG = (CP(ItJ) + CP(K,J))/2.0 

C CALCULATION OF EFFECTIVE FORCE 

FORCE = PAREAACPAVG 

IF( v .EQ.l ) GC TO 1 5 C 
IF ( J .EQ.NANG) GO .TO 150 
FORCE = 2 .0*FCRCE 

C CALCULATION CF MOMENT ARM 
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1 5 C XP = ((PC + RI)/2.C)*CCCS(THETA(J) ) 

C CALCULATION CF MOMENT 4 E 0 U T Y-AXIS 

YMCN = FORCE*XM 

C SUMMATION CF TOTAL FORCE 

TCTFCP = TOTFCR + FORCE 

C SUMMATION CF TOTAL MOMENT 

1 Yf^CN = TYMCM + YMCM 

160 CONTINUE 
1 7 C CONTINUE 

C CALCULATION CF CENTER CF ACTION 
X E A P = TYMCM/ TOTFOR 

C CALCULATION CF NCRMALI ZEC INUUCEO FORCE: Y/T 

\RINV = 1.0/VP 

YT = 2.0*TCTFCR*VRINV**2/PI 

C DATA CLTPUT 

V\RIT£(6»260) C8L0CK ' 

URITE<6,220) TCTF0P»X6AR 
V» R I TE(6»240) YT, VRINV 

180 CC 220 I = ltNANG 

ANGLE = THETA (I)*18C. O/PI 
V>RITE(6,190) ANGLE 

19C FCPNAT( *0' ,2X , 'THETA = » ,F6 . 2 , 8X , • R • , 9X, ' CP ' . 9X * • L ' , 
X 9 X » * V ' ) 

CC 210 J = 1 » NR AD 

VR I TE ( 6 ,2 CC ) P(J),CP(J«I),U(JtII tV(Jtl) 

2CC FCPNATC • ,23X,F6.2,2X,F8.2,3X,F6.2,4X,F6.2) 

210 CCNTINUE 
22 C CONTINUE 

230 FCRN AT (' O' »2Xt' RESULTANT FORCE OF • , IX, F7.2, IX , 

X ' ACTS AT X = • ,F6.2 ) 

240 FGRMAT('0',2X,'Y/T = • ,F6.2 ,1 X , 'FCR UM/UJ = ',F6.3> 
kRnE(7,25C)<XP<I),YP(I),ZP(I),I = l,MF0INT) 

250 FORM AT ( 3F20 .4 J 

260 FCRNAT(' O' »2X»' EFFECTIVE BLOCKAGE CYLINDER' , 

X'C IAMETER = • ,F4.2 ) 

27 C STOP 

ENC 
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